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This  technical  note  describes  cohesive  sediment  erosion  by  the  entrain¬ 
ment  process,  which  can  be  important  to  the  open-water  dispersion  of  flu¬ 
idized  dredged  material.  Some  specific  predictive  algorifiuns  are  presented. 

Background 

Entrainment  occurs  when  a  turbulent  fluid  layer  incorporates  fluid  from 
an  adjacent  layer  into  itself.  Entramment  is  a  fluid  process  that  can  be  im¬ 
portant  to  cohesive  sediment  dispersion  in  the  following  dredging-related 
situations: 

•  During  the  convective  descent  phase  at  disposal. 

•  Just  after  disposal  at  fire  leading  edge  of  a  rapid-spreading  rmderflow. 

•  Following  disposal  if  a  residual  fluid  mud  layer  is  formed. 

•  After  fluidization  of  disposed  sediments  by  waves  or  other  conditions. 

•  In  navigation  channels  near  agitation  or  overflow  dredging. 

Situations  in  which  entrainment  occiurs  can  thus  be  grouped  into  two  cat¬ 
egories  depending  on  whether  entramment  occurs  at  horizontal  or  vertical 
interfaces.  The  first  two  situations  above  represent  the  vertical  category. 

All  other  dredged  material  cases  involve  stable,  approximately  horizontal 
interfaces.  This  distinction  determines  whether  density  differences  control 
the  respective  entrainment  and  whether  the  entrainment  is  considered  an 
erosion  process. 

Entrainment  involving  horizontal  density  layers  is  an  important  erosion 
process  that  can  mobilize  disposed  sediments  for  transport  offeite.  Entrain¬ 
ment  in  cases  of  vertical  interfaces  is  not  inhibited  by  density  differences. 
The  general  ctmcept  of  entraiiunent  can  be  applied  to  a  nun^r  of  related 
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geophysical  processes  and  has  been  used  in  models  such  as  the  Corps' 
STFATE  disposal  model  and  other  near-field  models. 

Erosion  by  entrainment  is  controlled  by  density  differences  between  the 
layers,  particle  settling,  and  viscosity  of  the  stispension  as  well  as  flow  con¬ 
ditions.  Entrainment  can  occur  when  cohesive  sediments  are  in  a  flxiid 
state,  sometimes  called  fluid  muds,  with  relatively  low  concentrations  rang¬ 
ing  ^m  about  5  to  300  g/L,  depending  on  sediment  properties.  At 
higher  concentrations,  sediment  properties  change  rapidly  and  other  ero¬ 
sion  processes  come  into  play,  such  as  those  described  in  Dredging  Re¬ 
search  Technical  Notes  DRP-1-07  CTeeter  1992a).  At  lower  concentrations, 
settling  properties  of  cohesive  sediments  will  not  establish  a  sharp  luto- 
cline — a  region  of  high  vertical  gradient  in  sediment  concentration. 

The  Dredging  Research  Program  (DRP)  Technical  Area  "Analysis  of 
Ehredged  Material  Placed  in  0^^  Water"  is  testing  entrainment  process 
descriptors,  and  enhancing  them  as  necessary  to  accurately  describe 
dredgki  material  entrainment.  This  technic^  note  presents  background  in¬ 
formation  on  the  entrainment  process  as  it  applies  to  dredged  sediments 
and  describes  intermediate  results  in  algorithm  development. 

Additional  Infomiation 

Contact  the  author,  Mr.  Allen  Teeier,  (601)  634-2820,  or  the  manager  of 
the  Dredging  Researdi  Program  (DRP),  Mr.  E.  Clark  McNair,  (601)  634- 
2070,  for  ad^tional  information. 

Introduction 

A  fluid  mud  layer  lying  at  the  bottom  of  a  water  body  can  be  defined 
as  having  three  regions  in  the  vertical  dimaision.  Vertical  layering  above 
a  given  dense  layer  suspension  is  shown  in  Figure  1  for  high-  and  low- 
current  cases,  llie  upper  region  is  called  the  mixed  layer  and  is  assiuned 
to  be  turbulent.  The  intermediate  layer  is  called  the  stable  layer,  lutocline, 
or  interfodal  layer;;  its  presence  dep^ds  on  the  entrainment  process.  The 
intermediate  layer  may  move  with  the  flow  and  be  at  least  partially  turbu¬ 
lent.  The  lower  layer  is  called  the  dense  layer.  It  is  stationary  or  slowly 
moving  as  a  laminar  flow.  Because  of  their  relatively  low  settling  veloci¬ 
ties  and  cohesive  properties,  fine-grained  sediments  can  remain  in  a  fluid 
state  and  susceptible  to  entrainment  for  long  periods.  Surface  waves  can 
fluidize  cohesive  sediments,  as  described  in  [hedging  Research  Technical 
Notes  DRP-1-14  (Teeter  1994). 

Net  entraiiunent  through  the  lutocline  will  depend  on  conditions  in 
both  the  mixed  layer  and  the  dense  layer  and  thus  will  generally  be  time 
and  space  dependent.  The  net  entrainment  (E)  across  the  fluid  mud 
dense  layer  interface  is  the  difference  between  the  upward  turbulent  (f ^) 
and  downward  settling  (F^)  vertical  flux  components.  Thus, 
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Figure  h  Vertical  structure  above  a  fluid  mud  layer 
for  low-  and  high-current  cases 


E  =  F,  -  F, 


for  F,  >  F, 


The  turbulent  component,  F^,  is  initiated  by  motions  generated  by  d\e 
mixed-layer  flow  which  penetrate  die  lutocline  layer,  and  is  the  flux  rate 
at  which  sediment  is  mb^  upward  horn  the  dense  layer.  (A  notation  list¬ 
ing  is  provided  as  Table  1.) 

One  of  the  main  controls  for  the  turbulent  component,  F^,  is  the  overall 
density  difference  between  the  mixed  and  dense  layers.  The  relative  buoy¬ 
ancy  change  across  die  interface  is  Af’  =  g(p„  -  p^/p„  for  dense  fluids  and 
^  =  S(^m  ■  Cfc)(P  ■  Pi)/PsPo  suspensions,  where  g  is  the  acceleration  of 
gravity.  Varial7les  p„,  P{,  and  p.  are  the  densities  of  the  fluid  mud, 

mixed  layer,  fluid,  and  solid  partied,  respectively,  and  C„  and  are  the 
solids  concentrations  just  below  and  above  the  interfadal  layer.  Therefore, 
for  suspensions,  the  density  difference  (in  kilograms  per  cubic  meter  or 
grams  per  liter)  between  the  mixed  and  dense  layers  is  normally  about 
0.61  times  the  sediment  concentration  difference  in  grams  per  liter. 

A  velocity  depicting  the  rate  of  interfadal  displacement,  Vf,  due  to  tur¬ 
bulent  entrainment  is  ^/(Pm  ~  Po)  dense  fluids.  For  suspensions,  Vf 
=  Ff/iC„  -  Cj^,  where  F,  represents  the  flux  of  solids.  The  turbulent  veloc¬ 
ity,  Vf,  can  be  scaled  by  a  characteristic  velodty,  for  example  the  mean 
ntixed  layer  velodty,  U,  to  form  a  nondimensio^  entraimnent.  An  inter¬ 
face  displacement  rate  due  to  settling  is  defined  as  Ws,  =  (WsCj,)/(C„  -  Cf) 
where  Ws  is  the  settling  rate  and  WsC^  is  the  settling  flux  just  above  the 
interface.  Then,  a  nondUmensional  entrainment  for  a  suspension  is  (Vf 
-  Wsj)/U,  for  example. 
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Table  1 

Notation  Definitions 


Notation  Definition 


Amplitude  of  an  oscillating  grid,  m  _ 


Solids  content  at  the  bottom  of  the  mixed  layer,  kg/ cu  m  or  g/L 


Solids  content  at  upper  edge  of  the  dense  layer,  kg/cu  m  g/L 


Mixed  layer  depth,  m 


Netentiaimnent,  kg/sq  m/sec _  _ 


Upward  turbulent  flux,  kg/sq  m/sec 


Downward  settling  flux,  kg/sq  m/sec 


Gravity,  m/sq  sec 


Interfadal  layer  thickness,  m 


Viscosity  criteria  constcints 


Turbulent  entrainment  constant 


Vertical  eddy  diffusivity,  sq  m/sec 


Depth-averaged  Kz  for  nonstratified  conditions,  sq  m/sec 


Reynolds  number  for  the  mixed  layer 


Gradient  Richardson  munber 


Richardson  number  scaled  by  U 


Richardson  number  scaled  by 


Mixed  layer  current  speed,  m/sec 


Mixed  layer  shear  velocity,  m/sec 


Upward  entrainment  velocity,  m/sec 


Settling  rate,  m/sec 


Interfacial  displacement  rate  due  to  settling,  m/sec 


Vertical  coordinate,  m 


Buwancy  difference 


Mud  viscosity.  Pa  ■  sec 


Viscosity  of  mixed  layer.  Pa  ■  sec 


Mixed  layer  density,  kg/cu  m 


Fluid  density,  kg/cu  m 


Density  of  mud,  kg/cu  m 


Particle  density,  kg/cu  m 


Shear  stress.  Pa 


Angular  frequency,  radians/ sec 
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The  remainder  of  this  techiucal  note  will  address  entrainment  as  an  ero¬ 
sion  process  of  cohesive  sediments.  Four  mechanisms  important  to  the  en¬ 
trainment  of  cohesive  suspension  layers  are  described  in  Ae  sections  that 
follow:  density  stratification,  viscosity,  settling,  and  layer  thickness. 

These  mechanisms  act  to  control  entraiiunent,  while  the  entraiiunent  pro¬ 
cess  is  the  net  effect  of  multiple  mechanisms. 

Density  Stratification  Effects 

Upward  turbulent  flux  of  material  results  from  instabilities  and  recoils 
at  the  upper  surfiice  of  the  lutocline  and  from  intermittent  turbulence  and 
internal  wave  motions  that  reach  the  dense  layer.  For  horizontal  layers, 
density  differences  control  how  mudt  energy  goes  into  mixing  in  the  ab¬ 
sence  of  other  mechanisms.  The  greater  the  density  difference  between 
tibe  layers  in  relation  to  the  inertial  forces,  the  lower  the  £.  In  simple 
cases,  the  turbulent  entrainm^t  component  can  be  described  using  an 
overall  Richardson  number,  Ri,  relating  buoyancy  and  inertial  forces.  The 
most  appropriate  scaling  for  Ri  numbers  is  fixe  mixed  layer  depth,  D,  and 
mean  vdodty,  li,  although  for  mixed  layers  without  mean-shear  (su^  as 
tanks  with  oscillating  grids),  the  shear  velocity,  U*,  is  used.  Then,  for  ordi¬ 
nary  flowing  and  purely  tir^ulent  mixed  layers. 


AbD 

U2 


and  Ri* 


AbD 


(2a,  b) 


Equation  2b  applies  to  experiments  where  artificially  generated  turbu¬ 
lence  is  used,  as  will  be  desaibed  later. 

When  settling  and  viscous  effects  are  absent,  fluid  mud  is  entrained  in 
the  same  manner  as  a  dense  fluid.  Hgure  2  shows  entrainment  data  ob¬ 
tained  by  Naiimousa  and  Fernando  (1987)  and  Srinivas  and  Mdita  (1989) 
using  similar  racetrack-type  flumes  wher^  flows  were  driven  by  disk 
pumps.  The  Narimousa  and  Fernando  (1987)  data  are  for  salt  layers, 
while  flie  Srinivas  and  Mehta  (1989)  data  are  for  kaolinite-day  suspension 
dense  layers.  Kaolinite  is  a  common  commercial  clay  mineral  fluit  is  not 
particularly  cohesive  and  whose  settling  rate  is  much  lower  than  most  nat¬ 
ural  cohesive  sediments. 

hx  flxe  case  of  the  kaolinite  suspension  data,  the  Ws,-  is  very  small;  there¬ 
fore,  {Vf  -  WSf)/U  is  equivaloxt  to  V^/U.  The  Srinivas  and  Mdxta  (1989) 
suspensions  had  concentrations  of  50  to  133  g/L  clay  mixed  with  tap 
water.  The  mixed  layer  mean  velocity,  U,  ranged  from  50  to  100  mm/sec, 
wdxile  Ws  in  flxe  dense  layer  was  prdMbly  on  flxe  order  of  0.002  mm/sec 
or  less.  The  stable  layer  flxickness,  h,  with  respect  to  the  depth,  D,  was 
about  the  same  for  both  the  Srinivas  and  Mehta  (1989)  and  Narimoxisa 
and  Fernando  (1987)  tests. 
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Figuie2.  Entniiunent  for  a  dense  fluid  and  fluid  mud 
as  a  function  of  Riu 


The  line  in  Figure  2  associated  with  tt\e  Narimousa  and  Fernando  (1987) 
salt  solution  data  is  a  -1  power  law  fit  by  the  original  author,  and  since 
no  settling  was  associated  with  the  salt  ^ta,  die  observed  nondimensional 
entrainment  is  equivalent  to  Vf/U.  The  line  associated  widi  die  Srinivas 
and  Mdita  (1989)  data  is  a  regression  fit  to  the  log-transformed  total  en¬ 
trainment  data,  which  also  shows  a  -1  power  dependence  of  (V,  -  Ws,)/U 
on 

As  shown  in  Figure  2,  (V^  -  Ws,)/ li  or  Vf/U  for  both  data  sets  had  sim¬ 
ilar  functional  relationships  with  at  least  over  most  of  the  range. 

A  few  of  the  high-Ki„  kaolinite  data  points  have  reduced  entrainment  val¬ 
ues,  and  it  was  reported  that,  at  of  about  25,  turbulent  entrainment  at 
the  interface  ceased  as  the  result  of  suppression  of  turbulence.  Overall, 
the  data  in  Figure  2  represent  c(mditions  of  entrainment  where  suspen¬ 
sions  were  unaffected  by  settling  or  suspension  viscosity.  The  lines  for 
salt  and  kaolinite  density-stratified  cases  shown  in  Figure  2  are  similar  be¬ 
cause  the  entrainment  process  was  controlled  primarily  by  density  differ¬ 
ences. 

In  die  Ri  range  normally  involving  fluid  muds,  various  general  func¬ 
tional  forms  have  been  proposed  on  theoretical  grounds  and  supported  by 
various  data  sets,  the  most  widely  used  being  power  laws  with  -3/2  or  -1 
slc^Tes.  Arguments  for  a  particul^  functional  relationship  are  based  on  ei- 
tiier  energy  consideratiems  or  eddy  impingement.  Work  is  performed  by 
die  mixed  layer  as  dense  material  is  entrained  vertically  and  the  potential 
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energy  of  the  system  is  changed.  Based  on  data  presented  in  Figure  2 
and  by  odiers,  it  is  proposed  that 


F,  =  K(C„  -  Q  URi-'  for  P„UD/h„  >  k2 


(3) 


where  K  is  a  constant,  11^  is  the  dynamic  viscosity  of  the  mud,  and  the  sec¬ 
ond  expression  is  a  viscosity  criterion  based  on  a  constant  k2,  which  is  de¬ 
scribed  in  die  section  "Mud  Viscosity  Effects." 

At  low  Ri  values,  the  lutocline  can  become  largely  turbulent,  as  de¬ 
picted  in  die  hig^-current  case  in  Figure  1.  Turbulent  flux,  F,,  is  generally 
time  dependent  and  is  a  function  of  vertical  turbulent  ditfusion  and  verti¬ 
cal  concentration  gradient  values,  as  described  by 


F,  =  Kz(z)f 


(4) 


where  Kz  is  the  vertical  eddy  ditfusivity  and  the  parendietical  z  indicates  a 
value  at  a  given  point  in  die  vertical.  These  conditions  must  normally  be  cal¬ 
culated  numerically  using  various  computer  algorithms.  For  example,  die 
vertical  diffusion  can  be  ^culated  according  to  Teeter  and  Pankow  (1989)  as 


Kz(z)  =  /(z)  Kz„  [1+33  W^(z)l'-®  (5) 


where/(z)  describes  the  vertical  distribution  of  Kz(z),  Kz^isa  depth-averaged 
vertical  diffusion  value  for  homogeneous  conditions,  and  Rig  is  the  local  gra¬ 
dient  Richardson  number  given  by 


-gdpf,/dz  (6) 

Ri»=  — - ~ 

^  Pj,(divaz)^ 

Wolansld  and  odiers  (1988)  used  a  similar  method  to  obtain  predictions 
that  were  in  good  agreement  widi  held  measurements. 


Mud  Viscosity  Effects 

The  turbulent  component,  F^,  can  be  modified  by  viscous  ejects.  If  the 
fluid  mud  layer  viscosity  exceeds  a  certain  level,  interfacial  instabilities 
and  turbulent  entrainment  flux  rapidly  diminish.  Entrainment 
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suppression  has  been  observed  to  occur  when  viscous  and  inertial  forces 
are  comparable  at  the  turbulent  scale  of  the  entraining  flow  (Campbell 
and  Turner  1986).  Wide,  open  channels,  for  instance,  have  length  scales 
for  vertical  turbulence  approximately  equal  to  the  flow  depth,  and  the  ve¬ 
locity  scale  is  equal  to  the  mean  flow  speed.  The  ituxed  layer  criterion  for 
entraiiunent  is  a  Reynolds  nrunber  for  that  layer. 


Re=  p^UD/q,  > 


(7) 


where  U  is  the  average  flow  speed  of  the  mixed  layer,  D  is  the  chaimel  depth, 
and  q^  and  q^^  are  the  d)mamic  viscosities  of  the  mixed  and  dense  layers,  re¬ 
spectively.  liius,  entramment  begins  to  be  suppressed  when  li  D  <  q^ . 

On  theoretical  grounds,  Campbell  and  Turner  (1986)  proposed  that  the 
viscosity  effect  can  be  described  by  two  values  (kl  and  k2)  of  the  parame¬ 
ter  k,  where  A:  =  p^  l/D/q^.  Complete  suppression  of  entrainment  oc¬ 
curred  at  kl,  and  partial  suppression  occurred  between  kl  and  k2,  with  a 
linear  transition  between.  C^pbell  and  Turner  (1986)  estimated  kl  =  7 
and  k2  =  70.  A  plot  of  their  data  is  shown  in  Figure  3.  Based  on  the 
present  reanalysis  of  these  data,  the  values  kl  =5  and  kl  =  85  seem  more 
appropriate. 


Figure  3.  Relative  cntiaiiunent  according  to  Campbeli  and  Turner  (1986)  [points]; 
for  their  suggested  kl  ■  7  and  k2  «  70  [dashed  line];  a  regression  fit  to  their  data 
[dotted  line];  and  for  the  proposed  kl «  5  and  k2  «  85  [solid  line] 
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Viscosity  constraints  can  be  placed  on  Ff  for  suspensions  as  follows: 


p„UD-klr\ 


m 


m 


K(C„  -  Cj,)  URi-^ ,  fcl  <  P„  UD/r]„  <  k2 


(8) 


f,  =  0,  p„UD/n„<kl  (9) 

The  viscosities  of  fluid  mud  vary  widely  with  its  density,  Pj^^,  or  con¬ 
centration,  While  adequate  information  edsts  to  propose  descriptors 
for  entraiiunent  suppression  by  viscosity  ejects  in  negatively  buoyant 
pliunes  and  dense  suspension  layers,  the  most  appropriate  viscosity  charac¬ 
terization  has  not  yet  been  confined.  Fluid  mud  is  a  shear-thinning  mate¬ 
rial,  usually  exhibits  thixotropy,  and  can  be  fluidized  by  pressure  fluctua¬ 
tions.  Over  narrow  ranges  of  ^ear  rates  and  shear  stresses,  the  viscosity 
of  a  given  fluid  mud  can  vary  by  orders  of  magnitude  {Dredging  Research 
Technical  Notes  DRP-02-04,  Teetar  1992).  Density-viscosity  values  have 
been  found  to  vary  widely  between  sites,  with  temperature,  and  as  the  re¬ 
sult  of  other  foctors. 

Figure  4  shows  a  plot  of  observed  viscosity  versus  shear  stress  for  fluid 


Figare4  Viscosity  of  three  Sabine  River,  Texasy^  sediment 
samples  over  a  range  of  shear  stresses 

mud  from  foe  Sabine  Entrance  Channel,  Texas,  at  three  densities.  Varia¬ 
tions  in  viscosity  for  a  single  material  are  large,  and  at  some  threshold 
shear  stress  value,  viscosity  drops  markedly.  For  foe  range  of  shear 
I  stresses  normally  encotmtered,  say  less  foan  1  Pa,  Table  2  provides 
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example  data  on  the  relationship  between  solids  content,  bulk  wet  density, 
and  viscosity  for  a  cohesive  material. 


24 

41 

81 

106 

163 


Table  2 

Example  Cohesive  Material  Characteristics 
p„,kg/cam  Ti„,Pa  sec 


1,040 _ 1 

1,050 _ 5 

1,075  10 


Water  has  a  viscosity,  "nQ,  of  about  0.001  Pa  •  sec.  Using  density  versus 
viscosity  data  such  as  the%  and  disposal  site  information,  die  nomogram 
shown  in  Figure  5  can  be  used  to  judge  the  extent  of  viscosity  eff^ts. 


RgoreS.  Ncmiogiam  showing  k-values  for  complete  entniiunent  suppression  Qc  a  5), 
oitiainment  redaction  (5  <  k  <  85),  unhindered  entrainment  (85  <  k)  as  functions  of 

UDandrim/pm 


Entrainment  in  cases  of  vertical  interfaces  is  not  inhibited  by  density  dif¬ 
ferences,  although  at  high  concentrations,  cohesive  sediment  viscosity  may 
inhibit  entrainment  in  diese  cases.  Sediment  characteristics  have  been 
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previoiisly  shown  to  entrainment  coefficients  and  rates.  For  in¬ 
stance,  negatively  buoyant  plume  entrairunent  rates  were  found  to  de¬ 
crease  rapidly  and  approach  zero  as  moishue  contents  approached  the  liq¬ 
uid  limits  for  ti\e  material  (Bowers  and  Goldenblatt  1978).  Those  previous 
results  can  now  be  interpreted  in  light  of  newer,  more  general  experimen¬ 
tal  and  theoretical  results  on  liquids  of  different  viscosity. 

Particle  Settling  Effects 

Fluid  mud  settling  counteracts  upward  entrairunent  as  described  earlier. 
The  settling  flux,  Fj,  depends  upon  properties  of  the  sediment  material, 
concentration  of  suspended  se<fiment,  and  other  conditions  in  the  mixed 
or  lutocline  layer.  The  settling  flux  can  be  calculated  for  any  vertical  loca¬ 
tion  as 


Fj.  =  Ws  (z)  Cj,  (z) 


(10) 


If  conditirms  are  steady,  F^  can  be  evaluated  at  the  bottom  of  die  mixed  layer 
for  use  in  Equation  1.  For  the  low-Ri  case  when  ti\e  lutocline  layer  is  turbu¬ 
lent,  Fg  may  vary  vertically. 

Settling  rates  of  cohesive  sediments  vary  widely.  Fluid  mud  settling  is  nor¬ 
mally  in  the  hindered  range;  that  is,  settling  rate  decreases  with  increa^ 
suspension  concentration.  Over  most  of  foe  concentration  rai\ge,  cohesive 
suspensions  setde  as  a  mass.  Setfong  rate  is,  foerefore,  a  property  of  foe 
overall  suspension  rather  foan  of  foe  individual  particles  it  contains.  Larger 
silt  particles  and  day  aggregates  or  floes  all  setfie  togefoer  because  interparti- 
de  interactions  prevent  separation  and  sorting  by  grain  characteristics. 

Wolanski  and  others  (1992)  showed  that  fluid  mud  settling  decreases 
drastically  if  foe  fluid  mud  layer  is  mixed.  Decreases  were  observed  in 
settiing  at  foe  interfiice  by  factors  of  3  to  10  in  foe  presence  of  small 
amounts  of  turbulence.  Mixing  or  shearing  disrupts  microchannel  fluid 
pathways  important  to  foe  upward  flow  accompanying  hindered  settling. 

Within  the  fluid  mud  concentration  range,  at  least  one  concentration 
value  can  be  found  at  which  settling  flux  is  a  local  maximum.  This  oc¬ 
curs  typically  at  5  g/L  or  greater  concentration  depending  on  foe  sedi¬ 
ment  Where  these  concentrations  occur,  interfaces  form  with  abrupt  verti¬ 
cal  concentration  and  density  change  where  internal  waves  can  form. 

While  it  is  not  known  exactly  what  role  these  discontinuities  play,  it 
seems  likely  foey  influence  ^  vertical  strqcture  of  lutoclines  and,  thus,  en¬ 
trainment  Settling  flux  maxima  are  also  sediment  dependent 

The  effects  of  settling  on  entrainment  have  been  detected  experimentally 
and  reported  by  Wolanski  (1972).  Salt-  and  kaolinite-stratified  entrainment 
data  from  osdUating-grid  experiments  are  shown  in  Figure  6a.  These  tests 
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Richardson  Ntmiber  Ri. 
(a) 


(b) 

Figure  6.  Original  observed  dense  fluid  and  fluid  mud  entrainment  rates  (a)  from 
Wolanski  (197Z)  and  (b)  turbulent  entrainment  components  estimated  by  adding  a 

constant  settling  term 
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used  a  chamber  where  turbulence  was  generated  by  an  oscillating  grid 
rather  than  by  a  shear  flow,  such  as  would  be  the  case  in  a  flume  or  open 
channel.  For  these  tests,  die  Richardson  number,  Ri^,  is  defined  based  on 
shear  velocity,  ,  instead  of  U,  and  =  coa ,  where  co  is  the  angular 
rate  of  oscillation  and  a  is  ihe  amplitude.  The  clay  suspensions  had  low 
viscosity,  and  suspension  concentrations  were  on  the  order  of  3  g/L.  As 
noted  earlier,  kaolinite  has  settling  rates  about  an  order  of  magnitude 
lower  than  many  natural  cohesive  sediments. 

The  observed  entrainment  for  the  kaolinite  suspension  shown  in  Figure 
6a  included  the  effects  of  settling  and  was  much  lower  than  for  die  salt  so¬ 
lution.  Furthermore,  it  was  reported  diat  £  approached  zero  at  Ri  =2 
as  a  result  of  settling,  so  that  settling  flux  was  on  die  same  order  as  turbu¬ 
lent  entrainment  flux,  F^,  at  diis  point.  The  lines  in  Figure  6a  associated 
with  the  data  points  are  the  original  author's  fit  to  the  data.  The  data 
show  that  a  power  law  is  a  reasonable  approximation  to  the  relationship 
between  E  and  Ri^.  There  is  an  appred^le  difference  between 
(V^  ~  Ws.)  /U  versus  Ri^  relation^ps  for  salt  and  for  kaolinite  suspen¬ 
sions  obtained  in  the  same  experimental  apparatus  which  can  be  ex¬ 
plained  by  settling  effects. 

The  effect  of  settling  can  be  demonstrated  by  adding  an  estimate  of 
Ws./U  to  the  origi^  observed  entrainment  (V^  -  Ws,  )/LZ*.  Tests 
on  toounite  suspensions  were  performed  at  concentrations  of  only  about 
3  g/L  Under  diese  conditions  the  density  or  concentration  jump  across 
the  lutocline  is  small,  and  Ws.  approach^  Ws.  Settling  rate  was  not  mea¬ 
sured.  However,  a  Ws  value  of  0.1  mm/sec  is  representative  of  die  maxi¬ 
mum  settling  rate  for  kaolinite  (which  was  likely  to  have  occurred  at  this 
low  concentration)  and  was  assumed.  Figure  6b  shows  a  plot  of  how  the 
salt  and  kaolinite  suspension  data  lines  adjusted  to  yield  V  /U  .  Since  for 
salt  Ws  =  0,  die  V^/U^  line  for  the  salt  solution  is  identical  to  me  line  in 
Figure  6a.  The  assumed  constant  settling  rate  brought  the  V^/U^  line  for 
die  kaolinite  suspension  in  fair  agreement  with  die  salt  tests,  mdicating 
that  the  turbulent  entrainment,  was  equivalent  for  both  series  of  tests. 

Stablf  2nd  Dense-Layer  Thickness 

Entrainment  across  a  density  interfece  causes  a  lutocline  to  form  at  the 
top  of  the  dense  layer,  as  previously  described.  The  thickness  of  die  salt- 
stratified  stable  layer,  h,  h^  been  found  to  depend  on  the  depth,  D,  of  the 
mixed  layer,  and 


h=  0.06  D 


(11) 


Under  these  conditions,  die  lutocline  resembles  the  low-current  case 
^own  in  Hgure  1.  This  relationship  was  confirmed  for  the  moderate-Rt 
kaolinite  entrainmoit  shown  in  Hgure  2  (Srinivas  and  Mehta  1989).  How¬ 
ever,  at  lower  Ri  values,  Wolanski  and  odiers  (1988,  1992)  showed  that  the 
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lutocline  layer  thickness  during  fluid  mud  entrainment  is  a  balance  be¬ 
tween  the  shear  velocity,  settling  rate,  and  overall  density  difference 
across  the  lutocline,  and  that 


h  = 


4L/J 

WsAb 


(12) 


describes  the  thickness  of  the  lutocline  layer.  Under  fliese  conditions,  the 
lutocline  layer  resembles  fl\e  high-current  case  shown  in  Figure  1.  The  6b 
here  refers  to  the  difference  between  the  lutocline  layer  and  the  mixed  layer. 
Field  observations  show  that  flie  lutocline  layer  can  occupy  more  than  half 
the  water  coluxrm  depth,  O.  Lutoclines  1  to  5  m  thick  have  been  observed 
where  currents  were  on  the  order  of  1  m/sec. 

Many  dredged  material  layers  are  initially  of  limited  thickness  after  dis¬ 
posal  For  instance,  the  fluid  mud  tmderflow  at  a  Tylers  Beach,  Virginia, 
pipeline  discharge  was  predicted  to  be  only  about  0.03  D  (Appendix  A; 
Ihevenot,  Prickett,  and  Kraus  1992).  Fluid  mud  flucknesses  are  often  too 
small  to  be  reliably  observed  in  the  Held. 

When  the  da:ise  layer  is  thin  relative  to  the  mixed  layer,  say  <0.06  D, 
vertical  velocity  fluctuations  may  be  suppre^ed  by  the  rigid  layer  below 
the  dense  layer.  Penetration  of  the  dense  layer  by  turbulent  motions 
might  also  be  suppressed.  A  set  of  tests  were  carried  out  to  test  tiiis  hy¬ 
pothesis  since  no  information  on  this  effect  was  found  in  the  literature. 

To  test  the  effects  of  fliin  dense-layer  thickness  on  entrainment  rates,  a 
series  of  tests  were  performed  using  the  rtewly  developed  particle  erosion 
simulator  (PES)  device.  The  DRP  adopted  this  device  to  characterize  cohe¬ 
sive  sediment  erodibility.  The  heart  of  flie  P^  is  a  cylindrical  chamber 
(5  in.  (12  cm)  in  diameter;  5  in.  high)  that  contains  an  oscillating  grid. 

This  general  type  of  device  has  be^  used  extensively  in  entrainment  ex¬ 
periments,  alAough  the  PES  has  additional  features  to  handle  fine-grained 
sediments.  The  will  be  described  in  more  4etail  in  a  future  D^ging 
Research  Technical  Note.  The  PES  was  used  to  test  salt  brine  dense  layers 
(density  of  about  1.19  g/cu  cm)  that  had  limited  thicknesses.  Tap  water 
was  us^  as  the  mixed-layer  fluid. 

At  the  beginning  of  the  tests,  the  dense  layer  was  injected  imder  the 
mixed  layer  in  the  chamber  and  the  water  depth  of  flie  mixed  layer  ad¬ 
justed  to  5  in.  The  inituil  thickness  of  the  dense  layer  was  varied  in  ten 
tests  from  1.5  to  0.4  in.  (3.8  to  1  cm)  or  from  0.3  to  0.08  D.  During  indi¬ 
vidual  tests,  grid  oscillation  rates  were  increased  in  steps,  dense  layer 
flUcknesses  were  observed  using  dye,  and  samples  of  the  mixed  layer 
were  taken.  From  precise  determination  of  salinity  of  the  mixed  layer,  en¬ 
trainment  rates  were  calculated  for  each  test  step.  The  experiment^ 


14 


Technical  Note  ORP-1-15  (April  1994) 


results  were  nondimensionalized  by  plotting  entrainment  per  oscillation 
versus  Ri^. 

Results  shown  in  Figure  7  indicate  that  there  was  no  discernible  effect 
of  dense  layer  thickness  on  entraimnent  rate.  Based  on  diese  results, 
dredged  material  layers  of  limited  thickness  would  be  expected  to  entrain 
as  described  earlier. 
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Rgure  7.  Entrainment  for  various  dense  la3rer  thicknesses  (numbers  indicate 
thickness  relative  to  D  asapercent) 


Application  of  Proposed 
Entrainment  Algorithm 

Density  stratification,  settling,  and  viscosity  can  be  important  entrain¬ 
ment  medianisms  as  shown  earlier.  Algorithms  proposed  for  entrainment 
under  various  conditions  can  be  used  to  make  estimates  of  dispersion 
horn  die  vicinity  of  a  disposal  site.  Information  on  the  density  of  the 
fluid  mud  layer,  its  settling  and  viscous  properties,  and  ambient  flow  con¬ 
ditions  at  the  site  are  required.  Estimates  of  dredged  material  density 
following  disposal  could  be  made  by  a  separate  model  for  near-field  dilu¬ 
tion,  such  as  tite  STFATE  model  described  in  Dredging  Research  Technical 
Notes  DRP-l-()4  (Johnson,  Tallent,  and  Fong  1992)  and  DRP-1-10  (Nelson 
and  Johnson  1992)  or  other  models.  Determination  of  sediment  properties 
generally  requires  laboratory  analysis. 

Elaborations  cm  steady-state  solutions  to  Equations  3,  8,  and  9  are  neces¬ 
sary.  For  tite  steady-state  assiunpticm  (where  h  and  £  are  constant  in 
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time)  to  valid,  entrainment,  £  ,  must  be  constant  vertically  across  the 
lutodine  Equation  3  was  developed  for  data  for  which  and  £  were 
ccmstant  in  vertical.  Since  concentration  and  F  will  vary  vertically 
across  the  lutodine,  F^  must  also  vary  for  £  to  be  constant.  (Also  F^  must 
have  a  minimum  value  that  is  greater  than  the  greatest  F  in  the  luto- 
dine.)  In  the  usual  case  where  the  dense  layer  concentration  exceeds  that 
at  which  maximum  settling  flux  occurs,  lutodine  settling  flux  can  be 
estimated  as  the  maximum  settling  flux  for  titat  particul^  material.  Use 
ti\e  maximum  settling  flux  to  compare  to  a  calculated  F^.  If  the  settling 
flux,  F^  ,  is  even  as  large  as  0.25  F^ ,  entrainment  might  increase  with 
depth  tn  ti\e  lutodine  and  the  steady-state  assumption  may  not  be  valid. 

Lutodine  thickness,  h,  should  be  estimated  using  Equation  IZ  If  the  re- 
sidt  is  much  greater  than  0.06  O,  Eqtiations  3,  8,  and  9  should  be  used 
with  cautidL  Numerical  solution  methods  for  Equations  4  and  5  might 
be  required. 

When  the  greatest  level  of  accuracy  is  needed  or  other  circumstances 
dictate,  tiie  proposed  algorithm  can  be  used  as  part  of  a  numerical  solu¬ 
tion  sdieme  that  indudes  feedbacks  between  conditions  of  tiie  mixed,  sta¬ 
ble,  and  dense  layers.  Such  procedures  are  being  ir^nrporated  into  a  DRP 
personal  computer  program  now  tmder  developmeii 

Thougjh  only  one  application  has  been  made  tims  fax,  results  appeared 
to  overpredict  entrainment  son^what  (Thevenot,  Psrickeit,  and  Kraus 
1992).  The  proposed  entrainment  algorithms  wfll  continue  to  be  tested 
against  existing  field  information  and  any  new  monitoring  data  on  dense 
layer  behavior  in  the  field  that  becomes  avjulable. 

For  the  case  of  negatively  buoyant  plumes  from  continuous  or  instanta¬ 
neous  dredged  mata^  discharges,  entraiiunent  is  not  inhibited  by 
density  diffoenoes,  although  cohesive  sediment  viscosity  n.iay  inhibit 
entrainment  Descriptors  of  viscous  mechanisms  that  might  suppress  en¬ 
trainment  could  be  incorporated  into  numerical  model  formulations. 
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